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PREFACE 

This document reflects a complete revision of the orbiter powered flight 
guidance scheme. A unified approach to powered flight guidance has been taken to 
accommodate all phases of exo- atmospheric orbiter pow e re d flight, from ascent 
through deorfait. The guidance scheme has been changed from the previous modified 
version of the Lambert Aim Point Maneuver Mode used in Apollo to one that employs 
Linear Tangent Guidan e concepts. As such it supersedes the previous document. 
"Powered Flight Guidance". GN&C Equation Document No. 11. Rev. 2. April 1972. 
In addition, this document replaces the previous ascent phase equation document 
titled "Powered Ascent Guidance". GN&C Equation Document No. 1-71, January 
1971. 
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FOREWORD 


This document is one of a aeries of candidates for inclusion in a Mure 
revision of JSC-04217 "Space Sbuttle Guidance. Navigation and Control Design 
Equations". The enclosed has been prepared under NAS9- 10268, Task No. IS- A, 
"GNiC Plight Equations Specification Support", and applies to functions 1 through 
4 of the Orbit Insertion Guidance Module (OG1) and to functions 2 through 6 of the 
Poorer Flight Guidance Module (OG2) as defined in JSC-03690. Rev. D. "Space 
Shuttle Orbiter Guidance, Navigation and Control Software Functional Require- 
ments", dated anuary 1973. 




Gerald M. Levine 

Division Leader, Guidance Analysis 
NASA Programs Department 
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NOMENCLATURE 


a L 


Specific force limit during SSXIE maneuvers 


*T.t 


Assumed thrust acceleration for phase 



Constants used to detine tareet (entr> interface) 

velocity constraint, v * C. * C, 

v I . n 


C 


u 


Intermediate variable used to determine sensitivitv 

of v M to v 
—a t 


1 SS\1E 

! o\is 

'no 


Vacuum thrust of sineir SSME, OMS, or RC> engiu* . res pee u 


T.i 


Total assumed thrust for i T>l phase 


Gravity vector 


If Thrust acceleration integral, 

t 

/ COtl 

H- I '* 1 element of H. f ffj/r ») t“ dt 

l £Q, i-1 

Unit thrust vector 


t 

i co 2 

/ (f T /n ,)t «i: 

0 



Unit vectors relative to desired trajectory plane: 

is radial along r^. £ v is normal to desired trajectory 
plane, and ji^ is in the dov\ nrange direction. r J. x y -I * * 


— N 


I nit vector in the direction of the oroital angular 
momentum vector normal to the transfer plane (see Met. 5) 


-rt 


Unit vector in direction of £ ( 



i. j 


J 

J . 


K i 


K 

i 

k 

L 


max 


m 


Clement subscript variables 

/go 

Thrust acceleration integral, / (f^./m) t dt 

t 0 

fgO, i 

i 1 * 1 element of J. J (f T /m) t dt 

t go, i-1 

SSME throttle setting. i th phase (1.0 * 100%) 

Maximum throttle setting of SSME 

Subscript variable referring to current phase number 

} go 

Thrust acceleration integral, J (f_/m) dt 
t 0 
/g°. i 

ith element of L, J (f^/m) dt 

^o.i-l 

Current e. .imated vehicle mass 


m 0,i 

I *SSME 

^OMS 

^RCS 


Mass at beginning of i** 1 phase 

Mass flow rate of single SSME, OMS. or RCS engine, respectively 

Total mass flow rate, i th phase 
’.amber of thrust phases 


rev 


P 

P i 

Q 

Q,- 


Number of revolutions of coast (see Ref. 5) 

l go 


. / / (f T /n) s 2 ds dt 


Thrust acceleration ircegral, j j 

. 0 0 
fg °> 1 t 

i** 1 element of P, J f (f T /m) s ds dt 

^o. i-1 *go,i-l t 

Thrust acceleration integral, f / (t»/ m) s ds dt 


/go. i 

ith element of Q, J 




j (f T /m) s ds dt 
l go,i-l ‘go, i-1 
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£ 

— bias 


-cl 


— c2 


Id 


I go 


r 


coxy 


r 

goz 


r 

— grav 




I ref 


It 


Itef 

—thrust 


Vehicle position vector 

A position bias to account for effects of a rotating thrust v. v tor 
bias ~ — go " —thrust ^ 

Vehicle position vector at beginning of gravity computation 
coast segment 

Vehicle position vector at end of gravity computation coast 
segment 

Desired terminal (cutoff) position 

Desired radius magnitude at terminal (cutoff) position 

Position-to-be-gained including bias <£ thrust reflects true 
pos it ion*- to-be- gai ned ) 

Projection of r on plane defined by i and i 
—go r -x y 

Component of along (downrange) 

Second integral of central force field gravitational acceleration over 
thrusting maneuver 

Predicted terminal (cutoff) position 

Position on reference trajectory at tirre t ref 

Target position in inertial coordinates 

Entry interface target in Earth fixed coordinates 

Second integral of thrust acceleration vector over thrusting 
maneuver 
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engoff 


t t 

g° 

Thrust acceleration integral, J / (f„ /m 1 ds dt 

0 0 1 

Switch indicating whether engine-off command has 
. . , < - 0 , command not issued 

been issued . ( =1 command issued 


guess 


Switch set by Lambert routine indicating type of solution (see Ref. 5) 

r i 

ith element of S, J J (f^/m) ds dt 

^o.i-l ‘go.i-l 


’mode 


Maneuver mode 


1 Ascent, standard 

2 Ascent, reference trajectory 

3 Ascent, Lambert 

4 Ascent, once-around abort 

5 Ascent, return- to-launch- site abort 

6 On-orbit, external delta-v 

7 On-orbit, Lambert 

8 On-orbit, deorbit 


J - 1 first active guidance call 
P* ss l |=0 not first active guidance call 

s Gravity perturbation switch ( > 0 for non-Keplerian model) 


|=0, constant tnrust during phase i 
phase, i | a ^ constant acceleration during phase i 


Prethrust switch 


1=1, prethrust 
j = 0, active guidance 


Switch indicating whether the initial and target position 
vectors are to be projected into the plane defined by 
(see Ref. j) 

Lambert solution type switch (see Ref. 5) 


SSME , i Number of SSME, OMS, RCS engines, respectively, assumed 
OMS,i thrusting during i th phase 


RCS, i 


IX 



—bias 


Time associated with £ , v 

Estimated bum time remaining in phase i 

Coast time between last and next-to- last phase for tank separation 

Nominal cut-off time for reference trajcctor\ mode 

Time-to-go until end of maneuver 

t of previous call 

Time-to-go until end of ith phase 

Ignition time of first phase 

t of previous guidance cycle 

Time at target (point where terminal constraint defined) 

Iteration variable determined in Lambert Routine (see Ref. 5) 
Vehicle velocity vector 

A velocity bias tc account for effects of a rotating thrust 

vector (v, . = v - v , . ) 

-bias -go -thrust 

Vehicle velocity vector at beginning of gravity computation 
coast segment 

Vehicle velocity vector at end of gravity computation roast 
segment 

Desired velocity vector at terminal (cutoff) position 
Pern red velocity magnitude at terminal (cutofn position 
Value of v Q resulting from perturbed time of flight 


x 



ex, i 


-go 


goz 


V' 

-go 


V 

-grav 


-ref 


It 


It 


Ith 


v ' 

I th 


I thrust 


I tv 


V • 

- tv 




Effective exhaust gas velocity for phase i 

Velocity-to-be- gained including bias. thrust reflec; - 

true velocity-to-be-gained) 

Component of v along i (downrange) 

Recomputed to satisfy terminal constraints 

First integral of central force field gravity acceleration over 
thrusting maneuver 

Velocity on reference trajectory at time t re ^ 

Velocity vector at target 

Value of v t resulting from perturbed time of flight 

Projection of v ^ on horizontal plane 

Projection of o? horizontal plane 

First integral whrust acceleration vector over thrusting 
maneuver 

Projection of v ^ on vertical 
Projection of y' t on vertical 

Desired inertial flight path angle at terminal (cutoff) position 
Position offset used in gravity computation 


Ar z Downrange component of terminal position error (used during 

reference trajectory mode) 

At Guidance cycle time step 


Value of t used to define time to issue engine cutoff 
command and terminate active steering computations 
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At g° 

At to 




Av 

-go 


sensed 


Change in used in throttling computations 
Time interval before t^ to start active guidance calls 
Impulse velocity increment used in gravity omputations 
Change in v^ o 

Total velocity change accumulated on accelerometers since 
last reading 


AVQ^ig Total characteristic velocity to be imparted during the last 
ascent phase by the OBIS engines 

6t Perturbation in coast time used in deorbit re" :ired velocity 

computations to determine sensitivity to coast time 


c cone 

*vgo 

X 

X 

P 

a 

T i 


Sine of half cone angle of elusion (see Ref. 5) 

Value of lAy^l defining prethrust convergence limit 

Unit vector in direction of v ^ 

-go 

Time derivative of unit vector coincident with X but rotating 
wit. tesired thrust vector turning rate, u.*f 


Pam ping factor used in determining the change in Av^ q 


Scaling factor in required velocity computations 
Ratio of mass to mass flow rate for phase 



Time rate of change of 0 


Desired thrust .ector turning rate 
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INTRODUCTION 


1- 1 Objectives and Requirements 

The primary objective of powered flight guidance is to issue proper steering 
and if necessary throttle commands during the thrusting portions of a mission such 
that the desired objectives of the maneuver are satisfied in a reasonably efficient 
manner. In addition, the navigated vehicle state vector must be maintained through 
the maneuver. These are obvious requirements of the powered flight guidance pro- 
gram. Additional objectives or requirements, which influence the overall design, 
are listed below: 

(1) The guidancv program should be applicable to all powered 
maneuvers. This objective assumes that a single guidance 
routine is less expensive to design, code, and maintain in 

a flight computer than several smaller routines. This guid- 
ance scheme, referred to as Unified Powered Flight Guid- 
ance (UPFGh can be used for all exo-atmospheric orbiter 
thrusting maneuvers, from ascent through decrbit. Future 
versions of this scheme should probably include the atmo - 
spheric solid rocket boost phase for completness and to 
minimize any difficulties in transition from one guidance 
phase to the next. A complete description of the many types 
of maneuvers possible with this UPFG routine is included 
in the following section. 

(2) The guidance program should be simple and flexible. \ny 
guidance program which must handle the mam types of ma- 
neuvers required of the space shuttle orbiter. from the 
relatively high accelerations during ascent to the very low 
accelerations experienced with a single Orbital Maneuver. nc 
System (OMS) engine, can probably not be classified as 
truly simple. However, by properly structuring the routine, 
the various requirements can Lc handled efficiently and the 
impact of new requirements can be minimized. This UPFG 
routine is structured such that communality of basic com- 
putations is maintained with specialized operations per- 
formed according to the mission phase and desired maneuver 
objectives. In addition, to minimize program size, maxi- 
mum use is made of subroutines required for other G^N 
functions. 
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(3) The premaneir c 1 prediction of th ; maneuver Av require - 
ment should be reasonably accurate . In order to assist the 
premaneuver targeting process and make mission critical 
decisions, the powered flight guidance scheme must accur- 
ately predict the maneuver £x\ , especially for Long, low ac- 
celeration orbital maneuvers. This requirement can have a 
very significant effect upon the guidance algorithm design. 
Many guidance schemes are based upon approximations which 
become more accv! vt • as the man ::r. rt r progresses. Al- 
though these schem * may satis far - ' ily calculate steering 
commands during the course of the maneuver, they may not 
adequately predict the maneuver Av. This unifu c guidance 
scheme is based upon modifications to the Linear Tangent 
Guidance (LTG) concept which improve overall accuracy for 
both premaneuver targeting and guidance. This improvement 
in the LTG concept was also required for one of th** ascent 
guidance modes so that engine throttle changes, which are 
based upon an accurate prediction of the terminal (cutoff) 
state of the vehicle, could be properly calculated. 

(4) The guidance algorithm should satisfy primary maneuv er ob- 
jectives for nominal and perturbed conditions. When practi- 
cal, the guidance scheme should close the guidance loop 
around the true maneuver constraints, rather than a set of 
artificial constraints based upon nominal conditions. For 
example, during many on-orbit maneuvers the true con- 
straints lie on the coasting trajectory subsequent to the ma- 
neuver. During the terminal phase rendezvous maneuver 
the objective is to intercept the target vehicle. During the 
deorbit maneuver the true maneuver constraints are at entry 
interface, where a pre-determined relationship betw een entry 
range, velocity, and flight path angle must be satisfied. It 

is possible to determine a set of artificial constraints, de- 
fined at thrust cutoff, which will satisfy the maneuver objec- 
tives. However, the process of accurately caloulutii .* these 
artificial constraints can significantly complicate the target- 
ing process. In addition, the artificial constraints would 
necessarily have to be based upon a nominal coast trajectory 
alter thrust cutoff. \ny perturbations during the maneuver 
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could result in a Loss of performance (i. e. the true constraints are 
only partially satisfied) or an increase in total maneuver A v, 
caused by trying to force the solution onto the nominal coast 
trajectory at thrust cutoff. Furthermore, if the guidance 
system satisfies true maneuver constraint^ then much of the 
analysis necessary to determine guidance software perfor- 
mance under perturbed conditions can be eliminated. 

(5) 1 he guidance algorithm should produce near fuel optimal 

maneuvers for nominal and per rbed conditions. The Linear 
Tangent Guidance (LTG) Equations which form the basis for 
the UPFG scheme are based upon classical optimisation theory 
and appear to give excellent performance. However, it should 
be noted that even a truly optimal scheme will use excessive 
fuel if the ignition time is poorly chosen or the maneuver 
constraints result in an unnecessarily over-constrained ma- 
neuver which could be more efficiently performed by using 
more than one maneuver. 

1. 2 Types of Maneuvers 

From a guidance viewpoint, the orbiter maneuvers can be conveniently sep- 
arated into two classes. The first class, ascent, is characterized by a constraint 
on the vehicle altitude at thrust termination. The second class of maneuvers, on- 
orbit, does not require any constraint on the position at thrust termination. On-orbit 
maneuvers are typically intended to place the vehicle onto any coasting trajectory 
which satisfies a rendezvous intercept constraint or deorbi entry-interface condition. 
In these cases any type of thrust cutoff position constraint would probably increase 
the maneuver Av unnecessarily. 

Each of these classes of maneuvers can be further subdivided into individual 
maneuver modes with a particular set of objectives and a specified set of constraints 
which meet these objectives. These individual modes are listed below: 

(1) Standard Ascent Maneuver 

This mode is intended for use during the ascent phase of 
most miss ms. The insertion conditions are specified 
preflight, and defined by a desired terminal (cutoff) alti- 
tude, velocity, flight path angle, and orbital plane. Thus 
all components of tne terminal state are specified except 
the downrange component of position. 
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(5) Ascent to Coast Reference Trajectory 

This mode is intended for use on Mission 3B, a time 
critical mission involving ascent, rendezvous, satellite 
retrieval, and deorbit with return to the launch site one 
revolution after liftoff. To satisfy the objectives of this 
mission in the presence of perturbations and small launch 
delays, the main engine throttles are used to improve con- 
trol over the ascent and subsequent coasting trajectory. 
Introducing the throttle command as one additional degree 
of freedom in the guidance algorithm makes it possible 
to insert the orbiter onto a coast reference trajectory 
at cutoff. This coast reference trajectory can he deter- 
mined preflight to provide proper closing velocity with 
the satellite for the rendezvous braking maneuvers. 

f 3) Lambert Ascent Maneuver 

This mode can be used to insert the orbiter at a specified 
altitude onto a coasting trajectory which intercepts a 
specified position (target) at a specified time. It is 
similar to the coast reference trajectory mode, except 
that throttle commands for trajectory control are not 
used. Therefore, although the resulting coasting trajec- 
tory doej intercept the target, it does not guarantee 
the pr^ er closing velocity for sucessful braking. This 
mode car be used during the latter phases of ascent to a 
coast reference trajectory, when g -limiting considera- 
tions may override throttling. Since most per- 
turbs 4 Ions occur well prior to g- limiting, the coast 
reference trajectory mode is designed to complete com- 
pensatic n f .r perturbations prior to g -limiting. 

During k limiting, the Lambert ascent mode will main- 
v*' * «t intercept trajectory with only small deviations 
in terminal rendezvous closing velocity. It should be 
noted that no additional code is required to support this 
guidance mode, since it evolves as a natural result of 
two other modes, the on-orbit Lambert maneuver and 
the standard ascent. 
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(4) Pace -Around Abort 

This mode is intended for use during the latter portion 
of ascent in the event of an engine failure. The resulting 
coasting trajectory insures that proper entry interface 
condition 0 (range, velocity, and flight path angle) are 
achieved for sucessful reentry. As with other guidance 
modes, the insertion conditions necessary to satisfy 
true maneuver constraints, such as at entry interface, 
are recomputed every guidance cycle. Very little ad- 
ditional code is required to support this mode, since it 
evolves from a combination of the standard ascent and 
deorbit guidance modes. 

(5) Return-to-Launch-Site Auort 

(Studies are currently in process to determine desired 
thrust termination conditions). 

(6) External Delta- V Maneuver 

This maneuver is designed to guide the vehicle through 
a constant attitude maneuver which achieves a specified 
velocity change. This mode is used for small on-orbit 
maneuvers, such as rendezvous phasing maneuvers. It 
is similar to the Apollo External Delta-V maneuver 
mode. 

(7) Lambert Qn-orbit Maneuver 

This mode is designed to insert the vehicle onto a 
coasting trajectory which intercepts a specified position 
(target) at a specified time. This is typically a constant 
attitude maneuver and is intended primarily for rendezvous 
terminal phase and automatic braking . It is similar to 
the Lambert Ascent maneuver, except that no constraint 
is placed upon the vehicle altitude at thrust termination. 
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(8) Deorbit Maneuver 

This mode is designed to place the vehicle onto a coasting 
trajectory which satisfies entry interface conditions. These 
entry interface conditions are assumed to be defined by a 
prescribed (possibly functional) relationship between entry 
range* velocity* and flight path angle. This mode is similar 
to the once-around «:bort maneuver, except that no con- 
straint is placed upon the vehicle altitude at thrust termination. 

The functional flow of the program will be described in the next section to 
aid in understanding the general method employed in this guidance scheme. Follow- 
ing this* the particular equations and flow charts relating to the various modes are 
presented in greater detail. 
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2. FUNCTIONAL FUOW DIAGRAM 

JFFG is designed to be called by a Servicer Routine at appropriate 
times during a powered maneuver. The first call can be made at any time prior to 
active guidance calls. This first, or prethrust, call is different from the active 
guidance calls in that certain required parameters are initialized and the compati- 
bility of targeting information is checked and revised if necessary. The first call 
for active guidance is made at a specified time interval before ignition with following 
calls made at periodic intervals through-out the maneuver. 

During the prethrust call, an iteration takes place which recycles the rou- 
tine, without advancing the state vector, until convergence of the required velocity- 
to-be-gained takes place. During active guidance, however, a single pass through 
the routine is made each guidance c>*cle call, incorporating sensed velocity changes 
and updating accordingly. 

This section will present the general functional flow of computations per- 
formed l . UPFG. No attempt will be made here to define the actual equations used 
or to differentiate between the various modes of operation (other than to point out 
that several functional areas are not required by the External Delta- V mode). De- 
tailed equations and flow charts will be covered in later sections. 

The computational flow normally proceeds through nine distinct functional 
groupings or blocks of computations after each entry into UPFG. (see Figure 2-1) 
External Delta-V maneuvers, being less . implicated, by-pass two of these func- 
tional blocks. Each functional block of commutations results in the determination 
of the values of specific variables or commands used by the following functional 
blocks. The equations employed within the blocks may vary depending upon the 
various maneuver modes but the resulting output list from each block is the same. 

The first block of computations encountered is either initialization (Block 1) 
or update (Block 2) depending upon whether a prethrust or active guidance cycle 
call is being made. In Block 1 the initial values of required variables are set and 
the state is advanced on a coasting trajectory to the time when the first active 
guidance call will be made. In Block 2, after the accelerometers are read, the 
state is advanced based upon the sensed velocity change over the time elapsed since 
the last accelerometer reading. The sensed velocity change is also used to update 
certain other variables. Logic is included to exit the routine from this point if either 
(1) an active guidance call is being made prior to ignition or (2) an active guidance 
call is being made after an engine-off command has been issued. Prethrust flow 
does not use this exit path. 
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In the next block of computations (Block 3) the time-to-go remaining for 
the maneuver is computed. This computation involves solving for the phase elements 
of one of the thrust integrals. To avoid duplication these elements are saved to be 
summed .u the next block. 

Several scalar integrals relating to the thrust acceleration are computed in 
Block 4. These integrals are needed in following computations to predict the cut- 
off state and thrust direction. 

Block 5 incorporates an input thrust turning rate or determines one for ascent and 
abort modes. The turning rate is utilized to determine a unit thrust direction according 
to I.TG concepts. In addition, actual contributions to velocity and position due to 
incorporating the turning rate are determined and biases are computed which reflect 
the differences between actual contributions and the necessary velocity-to-be-gained 
and position- to -be- gained vectors. The biases are used in the routine to improve 
prediction. 

The steering block (Block 6» remains to be determined. It must, however, set up 
steering commands based upon the unit thrust direction . imputed previously, taking 
into ac ^unt autopilot requirements. This block of equations must diffet ntiate 
between prethrust and active guidance calls when making these computations. 

Block 7 and 8, which are skipped in the External Delta- V Mode, estimate 
the effects of gravity over the thrusting trajectory and incorporate these effects 
in the prediction of terminal conditions. Using these terminal conditions the 
velocitv-tc-be-gained is revised to force the predicted target conditions to satisfy the speci- 
fied target conditions. This revision is necessary to prepare for the next pass through 
the routine. 

Block 9, which determines and issues engine commands, is also to be de- 
termined at a later date when engine characteristics are better defined. It also 
must differentiate between prethrust and active guidance calls when formulating 
these commands. 

This completes one pass, prethrust or active guidance, through the r mtine. 

From this point, if a prethrust c?ll is being made and a change velocity-to-be- 
gained greater than some limiting value has been computed, the flow returns to 
the beginning of Block 3. If the velocity- to- be -gained change has converged to less 
than the limiting value, or if an active guidance call is being made, the flow exits 
to the Servicer Routine to await the next call. 
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ENTER 

l 


Prethrust 


Active Guidance 

i 


r _ _L 

Block l 

v 

Block 2 

• Initialise variables 


- Read accelerometers 

• Project state to t^- At tQ 


• Advance state 

• Update required variables 


(Prior to t. or after 
cutoff command) 
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3. 


INPUT AND OUTPUT VARIABLES 


The Unified Powered Flight Guidance program is not designed to run inde- 
pendently through-out the maneuver. It must be called periodically by the Servicer 
Routine, tust in prethrust and subsequently at the beginning of each active guidance 
cycle. Some variables require input values at each call and some variables are 
modified internally requiring only initial input values. Thus, the input list for 
UPFG can become complicated. Furthermore, certain variables optionally re- 
quire input values (initially and/or subsequently) depending upon the maneuver mode 
desired. In an effort to alleviate some confusion the input will be listed below in 
two catagories: input required on the prethrust call and input required each active 
guidance calL The first catagory will be further broken down into two groupings: 
those inputs required by all maneuver modes and those inputs optionally required 
by some modes and not by others. Variables that are modified internally only , and 
whose values affect future UPFG computations, will not be listed zs input. 

In addition to input which may be changed each call, a certain number of 
constants are required to be preset. These constants are not modified during 
program execution. They may be broken down into two catagories: universal con- 
stants and program constants. Universal constants are those constants not related 
to UPFG such as gravitational constants, earth radius, expected thrust levels, etc. 

Program constants are physical parameters which are unique to UPFG such as con- 
vergence criterion, iteration limits, etc. 

Sensors provide another source of input. The values of certain variables 
are updated by sensors on a continuous or semi-continuous basis. When a program 
uses the value of this variable it obtains the current value existing at the time the 
storage location for that variable is read. Sensors read and the variables they 
control are listed below. 

Many variables are modified and assigned new values during the prethrust 
and active guidance calls. With common storage locations they could all be termed 
"output". The list below, however, contains only those variables, known to affect 
other routines, that have been assigned new values by UPFG. Those that only 
effect subsequent calls to UPFG have been omitted. 

For on-orbit maneuvers (modes 6,7, and 8) the desired thrust vector turning 
rate u; £ is required as input to UPFG. It is anticipated that this quantity will be 
set to zero for External Delta-V and Lambert maneuvers, and computed in the de- 
orbit targeting program (Ref. 3) to minimize Av for the deorbit maneuver. If necessary, at 
a later date, an optimal turning rate for the Lambert maneuver could be determined. 
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M >| 


Input - Prethrust Call Only (All Modes) 


Symbol 


pre 


s 


mode 


Definition 
Pre thrust switch 


Maneuver mode 


i 1, prethrust 
| 0, active guidance 



Ascent, standard 
Ascent, reference trajectory 
\scent, Lambert 
Ascent, once- around abort 
Ascent, return- to-launch-site abort 
On-orbit, external & v 
On-orbit, Lambert 
On-orbit, deorbit 


t. 

ig 

* i 
i 

* 

s $SME. i 
s OMS.i 
s RCS, i 

3 phase, i 

V 

-go 


Number of thrust phases 
Ignition time, first phase 

State vector 


Number of engines (SSME, OMS, RCS) for i*h phase 

, 1 0, constant thrust 

Phase switch, l phase | constant acceleration 

Estimated velocity-to-be-gained -ectnr 

Mass at beginning of phase i 
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Input - Prethrust Call Only (Optional) 


Mode Required 


Symbol 


a 

B 

B 

a 

ES 

i 

D 

a 



B 

B 

1 

1 

1 

1 

B 

B 

B 





■ 

a 

B 


a 




■ 

a 

a 

B 

a 


1 



B 

B 

B 



i 



a 

B 

B 






a 

fl 

B 






B 

B 

B 

■ 





B 

B 
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OMS 



n 


rev 


S soln 



w f 


Definition 


Coast time between last and next to last phase 
Time remaining in i** 1 phase 
Velocity change to be imparted by OMS phase 
Specific force limit during SSME maneuvers 
Unit vector normal to desired tra ectory 
Desired cutoff position 
Desired radius magnitude of r ^ 

Desired velocity magnitude at cutoff 
Desired flight path angle at cutoff 
Desired velocity vector at cutoff 
Desired cutoff time 

Desired throttle setting for SSME for i** 1 phase 

Lambert unit normal to projection plane 
(see Ref. 5) 

Lambert number of revolutions (see 
Ref. 5) 

Lambert solution type switch (see 
Ref. 5) 

Target position - inertial coordinates 

Target (entry interface) - Earth fixed coordinates 

Time at target 

Constants in target (entry interface) velocity 
constraint, v v ' ^ ^2 v h 

Desired thrust vector turning rate 


Inputs for Mode 5 are to be determined. 
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Input Required - Xctive Guiri ice Call (All Modes) 


Symbol Definition 

^SSME k j q £ en g ines (SSME, OMS, RCS) to be considered on 

OMS, k ) for curren ^ phase (k*h phase) 

S RCS,k > 


Program Constants 
Symbol 


cone 


At 


to 


6 1 


At 


cutoff 


pert 


Definition 


Lambert required sine of half cone angle of exclusion 
Value of jAVg C | defining prethrust convergence limit 

Time interval before t.^ to start active guidance calls 

Offset in coast time used in deorbit requited velocity com- 
putations to determine sensitivity to coast time 

Value of t used to define time to issue engine cutoff 
command and terminate active steering computations 

Gravity perturbation switch 


Universal Constants 


Gravitational constant 


l SSME 

f OMS 

f RCS 


Fuf! thrust of single SSME, OM"\ Kt’S engine 


K Maximum throttle setting of SSMK 

max 


m 


SSME 


n^oyiS * Mass flow rate of single SSME, OMS, RCS engine at full thrus 


,fl RCS 


St 
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Sensors Required 


Sensed Variable Definition 

sensed Total velocity change accumulated oa accelerometers 

since last reading 

t Actual time when accelerometers are re*d (will be 

associated with state vector 

Output 

Output frc.a this program will be in the form of attitude and engine commands 
and state. These include: 

Definition 

Steering command (TBD) 

Ignition (engine-on) command (TBD) 

Throttle command 
Cutoff (engine-off) command (TBD) 

State vector 


Symbol 
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4. DESCRIPTION OF EQUATIONS 

The guidance scheme presented in this document evolved from a rather ex- 
tensive modification of the Linear Tangent Guidance (LTG) concept described in 
Ref* I. The original LTG concept was designed for an orbiter ascent maneuver in 
which the thrust cutoff altitude, orbital plane, velocity, and flight path angle were 
constrained to specified values determined prior to the flight. The current unified 
scheme retains that basic guidance mode, the standard ascent mode as described in 
the introduction* However, to adapt the LTG concept to the various shuttle maneu- 
vers, from ascent through deorbit, several important changes were made. 

First, the program was restructured to efficiently accomodate the various ma- 
neuver modes and to permit the calculation oi steering commands early in the guidance 
cycle, thus minimixing the computational lag. Second, the addition of a required 
velocity calculation toward the end of each guidance cycle was made. The basic 
standard ascent mode does not constrain either the downrange component of posi- 
tion or the time of thrust cutoff, and therefore the subsequent coasting trajectory 
can vary considerably depending upon ignition delays, atmospheric effects during 
boost, engine perturbations or failures, and any other factors which might effect 
the trajectory. By including a required velocity calculation every guidance cycle, 
based upon the predicted cutoff position, the true maneuver objectives can be satis- 
fied and the effects of perturbations minimized* The third important change to the 
LTG concept involved the elimination of any cutoff position constraint for on-orbit 
maneuvers. An alternate equation, which uses the input vehicle turning rate, simpli- 
fies premaneuver targeting and minimizes maneuver Av. The fourth change involved the 
design of a new scheme for the prediction of the effects of gravity on the powered trajec- 
tory. The original LTG scheme was not suitable for long, low thrust orbital maneuvers 
such as deoroit. A new scheme, based upon a conic coasting trajectory, appears to 
give jood performance for all msneuvc-s. Finally, to support the ascent to reference 
trajectory guidance mode, equations were added to account for a rotating thrust vec- 
tor. These additional equations improve the prediction of both the maneuver time and 
cutoff position. They are required for proper calculation of throttling commands 
during ascent, and they provide increased premaneuve. prediction accuracy in all 
modes. 


As described in Section 2, the Unified Powered Flight Guidance (UPFG ) 
Routine consists of nine blocks of in-line computations connected by the logic 
necessary to perform either the prethrust function or the active guidance function 
for the eight possible maneuver modes. A single entry and exit point are maintained 
for these computations. This section will describe the computations required for 
one pass, entry to exii, through UPFG. 
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Figure 5-1 shows the seque cing of the major computational block* Immedi- 
ately after entrv into L'PFG a test is made on the value of the switch, s . If 
J - pre 

Sp re ~ l, a prethrust call is being made; z 0 implies an active guidance call. 

For nrethrust calls, initialization in Block 1 is performed. For active guidance 
calls. Block 1 is bypassed and an update in Block 2 is performed. 


4.1 Initialization (Block 1) 


This blcck of computations provides the necessary one-time initialization of 
certain variables required to set-up the prethrust calculations. These variables 
include k. the maneuver phase counter. It should be noted that the first phase of 
each maneuver is denoted ”T\ Several switches are also set and the \alues of 
variables which may be defined by input, depending upon maneuver mode, are 
initialized. An initial value of £g ravr is required, since the calculation of r ^ ir. 
Block 5 occurs prior to the calculation of gravity effects in Block 7. To estimate 
— cra\ * assumed that the acceleration of gravity will be equal to its present 

value over the entire maneuver. Then, if a one second maneuver is assumed for 
simplicity. 


— gmv 


1 u I 

2 Ir* 


lja.ter, in Block 5, after the maneuver tin. o- to- go is computed, 
for the estimated maneuver time as follows. 


r is adjusted 

— grav J 


r (new) 
—grav 


t ^ (new! 
go 


r 


t (previous) 

g° * 


- 2f3V < t >reviOUs) 


where t (previous) has beer set to 1. On subsequent passes through Block 5 a 

previous t will have been set in the previous pass and the revision to r will 
r go —grav 

proceed normal tv. 


Finally, the effective engine paramete s for the array of engines to be used 
during each phase, i - 1 to n. are determined. These engine parameters include 
the total thrust (~ . , mass- flow- rate m . , the exhaust velocity v . , the 
initial acceleration a^. . , and the ratio of mass to mass- flow-rate r. for each 
phase of the maneuver. 


The last operation in Block l advances the <tate ti. r, v> from f hc nput 
values :•» the time when the '~irst actr e guidance .-all * Mt ne r arie at, t - t, - .it. 

This is acrot, p»i ^hed by calling the Precision State Rxtrapolati - Kc tine (Ref. 7). 
It should be noted that this sets up the prethrust operations as if Thrusting is ex- 
pected to commence at t. - At T(; instead of T . If is ^mail this is not 

expected to significantly affect tne results and it simplifies the program. An e.l- 
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tentative would be to extrapolate the state to t. in the prethrust call and each 
active guidance call prior to t^. It should be coted that this method of initializa- 
tion, although satisfactory for on-orbit maneuvers, is not applicable to ascent. 

In future revisions of this document a scheme which is satisfactory for both ascent 
and on-orbit maneuvers will be included. 


The computations now proceed to Block 3 since tne update (Block 2) is skip- 
ped during the prethrust call. 


4,2 Update (Block 2) 

Flock 2, which is skipped in the prethrust call, is the first block encountered 
during each active guidance call. This block acquires the velocity change sensed 
by the accelerometers and accounts for that velocity change by updating affected 
variables to the time when the accelerometers were read. The routine which reads 
the accelerometers and clock has not been defined but it is assumed that it 
will require some sort of initialization on the first call such that velocity changes 
on succeeding calls can be obtained by differencing from a previous reading. Thus, 
on the first call the velocity change is assumed to be zero. 

To allow for a variable guidance cycle time step, the time step is computed 
each cycle as follows. 

At - t - t 

prev 

where t is the current clock reading at the time die accelerometers are read and 

t is the previous time the accelerometers were read. On the first active 
prev 

guidance call At will be zero since *p rev was set in prethrust to - At^„ 

the expected time of the first active guidance call. The variable fp rev is set to t 
for use later. 

The velocity-to-be-gained is updated using the velocity change sensed by 
the accelerometers since the last guidance call as follows. 

y 5 V a Ay 

—go (new) -^o (previous) —sensed 

This gives excellent estimate of the velocity -to- be -gained for subsequent guidance 
calculations. Near the end of the guidance routine this value oi is adjusted 

slightly to account tor any changes in desired cutoff velocity and to prepai a for the 
next guidance cycle. 

It should be noted that on the first active guidance call (denoted by s^ssl " ** 
^Xs enS ed 13 assuine ^ *° ^ zero. Therefore, updating v^ Q is bypassed under these 
circumstances and 

= 0 
= 0 


sensed 


pass j 
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Next, the state, jr, v, is advanced to the time of the latest accelerometer 
reading by calling the Powered Flight Navigation Routine (Ref. 6). 

The remainder of the update block, which is bypassed when t £ t.^ , results 
in decrementing the t values for each phase, i, by At 

^o.i (new) ~ *go. i (previous) ~ 

during thrusting. Also decremented are the current mass, m, and current phase 
burn time 

m(new) = m(prevv'us) - rft k At 
T>,k (new) ’ S), k (previous) 


If t^ ^ ie^omes less than or equal to zero the end of that phase has been reached. 
Under ♦hose circumstances a change of phase is performed by setting the t . 
for the current phase to zero and then incrementing the phase number bv one. 

t. » 0 




go,k 
k(new) : 


k 'previous)^ I 


It should be noted at this point that mass, m, has not been adjusted to reflect step 
changes in this variable that may occur during a phase change due to such tilings 
as tank jettison or staging. It is expected that some routine external to UPFG 
will maintain current mass following phase changes for use by both Guidance and 
Control. 

During the r >ast period between the last and next-to-last phases of an ascent 
maneuver the coast time, t , is decremented by At. 


t c(new) " *c(previous) 


No changes to mass which would account for external tank jettison have been made 
here as this function ts assumed to be performed by a routine external to UPFG. 
Burn times remaining are also assumed to be constant during the roast time. 
Block 3 normally follows Block 2 during active guidance c Us unless t < t. or 
unless the :h ist-off command has been issued as indicated by s en g 0 ff = 1 • 

If either of these two conditions exists, all further computations are by- 
oassed and UFFG is exited. 
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4. 3 Time-To-Go 

The UPFG program is designed to provide guidance for all exo-atznospheric 
maneuvers, including ascent-to-orbit. The complexity of the time -to- go computa- 
tions is dependent upon the number of distinct thrust phases in the maneuver. All 
orbital maneuvers nominally have only one distinct thrust phase (OMS constant 
thrust) , while an ascent maneuver will have at least three distinct thrust phases 
(i.e. ( 1) SSME constant thrust prior to g-limiting , (2) SSME constant accelera- 
tion during g-limiting and ( 3) OMS constant thrust after tank separation). The 
equations that are described in this section are the equations required to compute 
the time-to-go until maneuver completion for a multi-phase ascent maneuver. 

The equations required for any single-phase maneuver are merely a subset of these 
equations; in particular, they are the same equations that will be solved during an 
ascent maneuver when actually in the final thrust phase of the maneuver ( OMS con- 
stant thrust phase) . 

lo compute the time-to-go it is first necessary to compute estimates of 
the current values of thrust magnitude f^. mass flow rate thrust accelera- 
tion ^ , effective exhaust velocity and mass to mass flow rate ratio, . 

where k is an index refering to the current thrust phase. Estimates of f^ and 
m^ are given by 

f Xk K k S SSME.k %5ME + S OMS. k f OMS * S RCS.k f RCS 

m k K k S SSME. k “*SSME + S OMS, k m OMS + s RCS. k “*RCS 

If £ 0, a^. ^ . v ftx ^ and r ^ are computed as follows 

a T. k f T.k^ m 
v ex, k f T.k^ “k 
T k = v ex.k/*T.k 

where m is the current estimated mass of the vehicle. If th^ = 0, the last three 
equations are bypassed and the previous estimates of a^, ^ and r ^ are 

used. 

If there are two or more thrust phases remaining (including the current 
phase) . the assumed bvrn times, . of each phase are used to compute the ve- 
locity change due to thrust, L. , that will be epplied during each phase i from the 
current phase to the n-2 phase (n denotes the number of phases in the maneuver). 

If phase i has constant thrust. 
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and if phase i has constant acceleration 


■i. 


b,i 


where a^ is the SSME acceleration limit. 


Since the velocity change to he applied in the n** 1 phase, L n , (after tank 
staging) is predefined, the velocity change to be applied in the n-1 phase can be 
determined by 


L i 



i- 


n-2 

L 

i=k 


L. - I. 

i 


where j - n-1 and v^ Q is the velocity- tr be-gained vector computed in Block 2. 
If there are onlv two phases remaining fk = n - l) t thrn 


t. 

J 


1 


v 

-go 


I * L 


n 


where again j n - 1. During the final phase of the maneuver L. is given by 


L i ' '-60 1 


where, in this case, j = n. 


Having determined L^, the burn time remaining in the phase can be 
computed. If phase j has constant thrust 

-I-./ v 


t . = T. (1 
Vj J 


J e\.j 


or if phase j has constant acceleration 


Vj 1 L i /a ^ * 

The time-to-go until the end of phase t . , for i 

go 1 1 

quired in Block 4 in order to evaluate the tiiru^; integrals L, 
These times are given by 


- k n is re- 

.1, S. Q, H. and P. 


t . = t . , + t tor i n 
go,i go. 1-1 b.i 


and where 


t . t * 0 for i * k 
go, i-l 


The time- tr .-go until the end of the n*h phase, which is the total 
until the end of the maneuver.is given by 


g° 


t go, n "go,n-l “b.n * 
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where t £ is the coast time during external task separation between phase n- 1 
and phase n. Although t c is not actually part of the maneuver time, it is added 
to the time-to-go until the end of phase n in order to maintain continuity of die 
steering commands. 

4.4 Integrals of Thrust 

The LTG guidance concept requires the evaluation of several thrust integrals. 
These integrals are defined as follows 


feo 

-d 

t 

I 

L * / (f/m) dt 

0 

( 

(f/ m) ds I dt 

rgo 

J = / (f/m) t dt 

0 

4 l 

t 

f 

. J . 

(f/m)s dsjdt 


H = / g °(f/m)t 2 dt P= (l / (f/m)s 2 ds dt 

0 J 0 l 0 J 


where (f/m) is the thrust acceleration, t is the time-to-go until the end of the 
maneuver, and t and s are variables of integration. For the space shuttle 
vehicle, it is assumed that either (f/m) is constant or f (thrust) and m (mass 
flow rate) are constant, therefore, these integrals can be integrated in closed 
form. 


The UPFG program is designed to accomodate multi-thrust- phase maneuvers 
as well as single -thrust- phase maneuvers. Therefore, the above integrals must 
be evaluated piece by piece since each thrust phase has a distinct thrust profile* 
This is accomplished by evaluating the thrust integrals separately for each phase 
and then summing them up. The thrust integrals for each phase i are defined 
as follows 


L i 


j. 


1 


H i 


*g°. i 

J (f T .i/ m i> ^ 


t 

go,i 

t 


S i = 

J 

j (f T i /m i )ds 

dt 

*go.i-l 


go. l-l 

*go,i-l 


*-go,i 


*go. i * 

f,‘ 

1 

J (f T ./m.)tdt 


/ 

[J (f T> ./m.)s dsj 

*go,i-l 


*80.1-1 
/ r 

*go,i-l 

/ 



A 

J “T.iW 2 dl 

*£ 0 . 1-1 


P i = 


go. i 


/ 


[/ (fT.i/mj). 


ds 


dt 
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where (f T . / m. ) defines the thrust profile for phase i , and t . is the time- 
to-go until the end of phase i (note: t 1 0 for i * k). 

The computation of the thrust integrals for each phase is performed in two 
steps. This is done in order to minimize computer memory requirements and 
computation time. In the first step, the equations vary depending upon the type of 
thrust phase, while for the second step, the equations are identical for both types 
of phases. Also, it should be noted that Hj is not explicitly comjwted. It can 
be shown using integration by parts that H, which is the sum of the H.' s . can be 
computed directly as a function of the time-to-go and the integrals J and Q. The 
evaluation of L. has already been described in Section 4. 3 because it is required 
in the computation of time-to-go. However, for the sake of completeness, it will 
be described again in this section. 


Step 1 in evaluating the thrust integrals is to comjwte L. and then to com- 
pute part of .1., S., Q., and P. . If phase i has constant thrust then 



P i * 




. + t 


go.i-1 




ex,i 


\> a i *3 *b,i 'go.i-l 


and if phase i has constant acceleration then 


L. = 

i 


l L Vi 


J i = 


2 L i Vi 


S. * 

1 


J- 


Q i 


F. 

l 


ViVi * Vi-' 1 


7T s. (t 2 . + 2 t . t . , + 3 t 2 . . ) 

6 l go.l go.i go, i-l go.i-1 
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In step 2, the remainder of J.. S., Q. and P. are computed as follows 



J i * L i 

V L Vi 
<*i + J Vi 
P i* H Vi 


where L, J, and U are the total thrust integrals from the current phase to the i- 1 10 phase. 
If i = n , the effects of the coast time, t , between phases n-1 and n are then added 
in as follows 


<*i ■ 

p i * V H ‘c 


Having evaluated the thrust integrals for each phase, the total thrust integrals 
are given by 


n 

L = E L. 

i=k 1 


n 

J = E Jj 

i=k 1 


n 

S = E S. 

i=k 1 


Q = 


n 

E 

i=k 




n 

P = E P. 
i=k 1 


H ’ l go 


It should be pointed out that H must be evaluated as the total thrust integral 
from the current phase to the i** 1 phase for each rliase because it is required in 
the computation of each . 



4. 5 Turning Rate 

The main results of this block of computations are the desired unit thrust 
direction, _if * and a vector, _X , which is associated with the thrust turning rate. 

The first operation in this block is to define X . a vector in the direction of 
Vg Q . the velocitv-to-be-gained. The unit vector, X, is the vector about which an 
expansion is later made to determine _i £. ■* From this point one of four different 
methods will be employed depending upon s m0( j e • 


-go * 


The first operation in this block is to define X . a vector in the direction of 
the velocity- to-be- gained. The unit vector, X , is the vector about which 
an expansion is later made to determine i_ f • From this point one of four different 
methods will be employed depending upon s mo< j e * 

Modes 1 through 5, ascent and abort, follow closely the LTG method of de- 
termining the rate X , by first using a desired burn-out position, r ^ , to estimate 


£go 


Igo 


r , - (r 4 vt + r ) 
— d — — go — grav 


The effect of gravity on this £^ q is given by £g rav which is estimated using £g rav 
from the previous call as follows. 


— gravfnew) 


2 

( go(new) \ 

*go( previous) I ~g rav \ 


previous) 


The projection of on the plane normal to the downrange direction, £g OX y 

is given by 


lz '- unit <Id x iy* 

— goxy —go z ’ — go ^ — z 


Using the integral, S, computed previously, the downrange component of can 

be modified by the LTG relationship 


goz 


is >»' 
X . i . 


— — z 


and a new r is thus found to be 
• go 


r = r f i + r 

-go —goxy goz -z —bias 


In this equation the effects of a rotating thrust vector are included by the term, 
r bias' wlu«'h was computed on the previous guidance cycle in Block 5. The rate , 
X , which corresponds to the velocity of the tip of a unit vector coincident with X 
but rotating with the desired unit thrust vector rotation rate is now obtained using 
♦he integrals, L, J, S, and Q: 

<r -S>) 

\ - KP ■ . >— 

- (Q - S J/L) 
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Note mat 


* * 

For modes 6-tt, a rotation rate, , is input and X is determined by 

X = <*2 £ unit [( X x r ) x X ] 

The predicted unit thrust direction at time, t, is given by 

2 f = unitfX - (J/L) X] 

It is recognizeu that the results of integrating a rotating thrust vector can 
be significantly different from r and v if the rotation angle is large. 

O o 

The angle o between X , a unit vector in the direction of Vg Q . and _if * the 
unit thrust direction, is given by 

d = cos’* n £ . X ) 


Since the linear tangent guidance equations are designed to align X and the thrust 
direction at the time J/L (J/L is approximately the midpoint of the maneuver), 
then 

0 = - d L/J 

Based upon 0 and ^ the first and second integrals of the thrust acceleration are 
given by 


—thrust 


—thrust 


jgo 


[ X cos(0 + dt) + sin(o + dt)] dt 

m - w 


r [ 4r 

0 0 


[X cos(d+ dt) +1^1 ?iti\ o x ot))di dt 


This can be simplified by asstiming that 

sin( d + ot) « d +0t 


and 


cos( O+dt)* 1 -( o ■** dt) 2 /2 


Using the thrust integrals computed in Block 4 , the actual (considering rotation) 
change in position and velocity due to thrust is computed as follows. 


thrust 


(L - | L <> 2 - J p + ^ H <i 2 ) X 

-(L d + J d ) unit (X ) 


"thrust = 

-(S ^ + Q 15 ) unit ( X ) 
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It may be noted that v thrust —thrust are reso ^ ve< ^ * nto components parallel 
and normal to X (parallel to X ). 

Biases to thrust cutoff velocity and position are computed by 


V - v “V 

—bias —go —thrust 

—bias —go —thrust 


4. 6 Steering Commands 

This block of computations is to be determined at later date. It will re- 
ceive the vector^ X and X and determine steering commands based upon the unit 
thrust direction as determined by an equation of the form 

j f * unit fX - (J/L) X] 

In order to issue steering commands to the autopilot, lead terms may be 
added as required. 

At the end of the prethrust call steering commands will take the form of a 
command to maneuver to the ignition attitude but during active guidance calls a 
turning rate may be implied. Logic will be included to incorporate this and to acti- 
vate and deactivate steering at the proper times. 
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4.7 Prediction of Gravity Effects (Block <# 

The solution of the LTG equations requires a prediction of both the first and 
second integrals of gravity over the ttuusting maneuver* The technique originally 
devised for ascent, and described in Ref. 1, was not appropriate for the long, low- 
thrust on-orbit maneuvers such as deorbit. Therefore a new technique has been 
devised which is applicaole to all maneuvers. This technique is based upon a 
coasting trajectory which is constructed such that it remains ’close’ to the powered 
trajectory thro- ghout the maneuver. The effects of gravity on the powered trajec- 
tory are then assumed to approximate the effects of gravity on the coasting trajec- 
tory. Thus the Kepler (Conic State Extrapolation) Routine can be used to deter- 
mine the required integrals of gravity. Figure 4-i illustrates this concept. 


Initial Coasting Trajectory 



Ignition 



Figure 4-1. Prediction of Gravity Effects 
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To construct this special coasting trajectory, assume for the moment that 
the maneuver takes place in field free space. The initial conditions ( t = 0) for the 
powered trajectory are defined by r and v * The a«st and second integrals of 
the thrust acceleration over the powered maneuver* i thrust afH * —thrust * are 
scribed in Section 4.5 . Therefore, at thrust cutoff he position and velocity on the 
powered trajectory are given by 

—cutoff ’ *go f —thrust 

—cutoff — * —thrust 

l 

where t^ is the maneuver time. The resulting trajectories in field free space 
are illustrated In Figure 4-2 9 where coasting trajectories simply appear as straight 
lines. 


‘Iti-Scl 1 



Figure 4-2. Prediction of Gravity Effects - Field Free Space 
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The initial state (t * 0, r. r) and cutoff state (t * t^. £ cuto jy. 1 cutoff * 
in field tree space are completely defined. A cubic equation can be used to model 
die state vector on die powered trajectory as a function of time. 

r p (t) = A + Bt + Ct 2 + Dt 3 
v p (t) = B + 2Ct + JDt f 

where r p (t) is die position. The velocity v p (t) is equal to d r p (t) / dt . The 
four vector coefficients. A. B. C, and B can be determined such that the boundary 
conditions on position and velocity at the initial and final times are satisifed. It 
should be noted that in actual practice it is not necessary to actually solve for 
these coefficients. They are merely used to aid in describing the concept. 

A coasting trajectory can now be constructed which remains 'dose* to die 
powered trajectory. The position on this coasting trajectory £ c (t) is defined by 
the linear equation 

£ c <t> = A* + B* t 

The velocity £ c <t) is constant and equal to B 1 . 

To determine A' and B' . the following integrals must be satisfied. 

(r c (tl - r pit) ] dt = 0 


t 

/' 


f go 


t c (t) - £ p (t) 1 ( tgo - t) dt 


0 


Since gravity is strictly a function of position, die first integral insures that the 
average position difference (or error) is zero. In addition, since errors in die 
initial position (and gravity) have more time to propagate and thus have more in- 
fluence on the total position er. or. the second integral weights the error as a func- 
tion of time. Using these integrals, the initial position A' and the initial velocity 
B' on the coasting trajectory can be easily computed. Using A* and B*. the initial 
conditions for the cc v-ting trajectory, _r cl and v reduce to the following simple 
form: 


r 

— c 


1 


v 

— c 


1 


— 10 —thrust " 30 -thrust *go 

v + 6. -thrust _1 

— 5 t 10 —thrust 

go 
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The Kepler Routine is usod to extrapolate these itucial conditions through the 
time t , thus obtaining r c2 and Then the effects of gravity on the coasting 

trajectory, which approximate the effects on the powered trajectory, are given by 

-grav — c“ —cl 
— grav " — c2 ~ —cl -cl *go 


4. 8 Velocity -to-be-Gained 


A more accurate prediction of the cut-off position is now obtained as follows. 

r s ^ v t r ■*> y 

— p — — go —grav —dirust 

For on-orbit Lambert and deorbit maneuvers ( s mo< j^ = 7,8) a desired 
cutoff position is given by 


id r i P 

For aborts and ascents (other than reference trajectory) the thrust cutoff 
altitude is constrained, therefore 

r , - r . unit (r ) 

— d d — p 


After determination of r ^ this block splits into three branches according to 
mode: standard ascent. ascent to reference trajectory and the remaining modes. 


4.8.1 Standard Ascent, - l 

The required velocity at burnout, v ^ , to satisfy the terminal constraints 
for this mode is a function of the inputs: J. v * r <j » - Since r d was alreadv 

used to determine , the remaining three variables ti vector and 2 scalars) 
are utilized in the following manner 


- unit (r^) 


±z = lx x —y 


Id = v d 


r i 1 

— X 

T 

sin> d 

i 


0 

— V 



i 


cos ^ . 

L — Z -J 


L d J 


Note that 2 X and i z define the radial and downrange directions with respect to iht 
cutoff state. The original LTG equations, described in Ref. 1. used the current 
vehicle position to define these directions. This resulted m continuous rotation 
the desired terminal velocity. The input normal to the transfer plane. * * s 
directed in the opposite sense to the current orbital angular velocity vector. 
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4.8.2 Ascent to Reference Trajectory jg m = 2) 

The objective of this mode is to intercept a coast reference trajectory for 
the single orbit rendezvous on Mission 3B. Throttling of the Space Shuttle Main 
Engines is used to increase the number of guidance control variables (i. e. , degrees 
of fr ee dom) by ooe, and thereby enable the guidance algorithm to control the down- 
range component of position at insertion. 


The following process is used to accomplish this. First, foe Conic State 
Extrapolation Routine (Ref. 4) is used to extrapolate foe reference trajectory state 
from foe input values (r to the predicted forest cutoff time t + t^ . (It 

should be noted that for proper operation of the guidance equations during the initial 


pass through the program, it is assu me d that the input reference state (r^, 
corresponds to foie nominal time of thrust cutoff determined from preflight simula- 
tion). The extrapolated state (r d . v^) can be used to determine foie error in the 
downrange component of position at cutoff (J z • (r ^ - r p) J. The other components 
of foie extrapolated state (r Vg ) are used to update v^ and r ^ on foie sub- 
sequent guidance cycle. 


Using this error in downrange position, it is necessary to determine a change 
in foie time-to-go which drives this error to zero. Based upon foie simplifying as- 
sumption of a flat earth and constant acceleration, it can be shown that foie change 
in relative position Ar £ resulting from changes in t^ 0 is given by 


d Ar 


d t 


go 


- i -v 
—71* _g° 
2 


Then the change in 
mined from 


go 


necessary to drive the position error to zero can be deter - 


At 


8 ° 


-2 



The new SSME throttle setting necessary to achieve this change in time-to- 
go is based upon the assumption that foie maneuver velocity-to-be- gained is fairly 
insensitive to changes in throttle setting. Thus 


K 


k 


*k *b.k 

Vk + ^go 


where L ^ is foe bum time remaining in the current maneuver phase. The throttle 
setting is limited to a maximum value K mJU£ . 
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On the last guidance cycle prior to the constant acceleration ( g- limited) 
phase, die maneuver is converted to die Lambert ascent mode = 3) and 

ascent is completed in that mode. This is necessary since throttling for trajec- 
tory control and g -limiting may be incompatible. However, use of an artifical 
g- limit sl-'ghtly less than die true limit, or a slight relaxation of the g- limit con- 
straint could make it possible (with minor changes) to maintain the reference trajec- 
tory throttling mode throughout die second maneuver phase prior to tank separation. 
The Lambert ascent mode, however, will insure intercept with the target satellite 
with only minor variation in closing velocity. 


4.8.3 Modes Requiring Lambert Solutions. = 3.4.7. 8 

The remaining modes (with the exception of s mo j e * 6, External Delta- V. 
which does not enter this block) require calls to the Conic Required Velocity Deter- 
mination Routine (Ref. 5). This routine solves a Lambert problem to determine 
the initial velocity required to satisfy certain terminal constraints. 

The straight Lambert modes, ascent and on-orbit = 3,7) require 

one call to the Conic Required Velocity Routine to determine, v^ . the desired 
velocity to coast from the end of die maneuver to the target, , in die given 
time interval. t| - (t + t ). 

Modes 4. S, 8, abort and deorbit require atmospheric reentry with the 
entry point at a constant location relative to die rotating earth. It is assumed 
that a constraint of die form mentioned earlier is required of die vertical, v^ . 
and horizontal, , components of entry interface velocity v ( . 

v tv =C l‘ C 2 v th 


For small entry angles (-1.0 to -In 7 degrees) and entry velocities around 26,000 
ft/ sec this constraint is nearly equivalent to expressing the entry angle as a linear 
function of the entry velocity. 

The only remaining variable that can be adjusted to cause the entry point to 
satisfy this terminal constraint is the coast time between cutoff and entry inter- 
face. t. * (t + t„ ). In other words, since t and t are fixed by previous com- 
t go go 

putations, t fc , the time at entry interface, must be adjusted. Using the currently 
* ssumed t^ and entry interface point relative to earth, r an inertial entry 
interface position, £ t , is determined by a call to the Earth Fixed to Inertial 
Routine. Then r ^ . r t , and the coast time (^ - t - t ) are input to the Conic 
Required Velocity Routine to obtain the desired velocity at cutoff, , and velocity 
at entry interface, v t . The entry interface velocity, v t , will not in general lie 
on the constraint defined above but at some other point as shown in Figure 4-3. 
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Horizontal velocity, v 



Figure 4-3. Entry Interface Constraint 
The components of Vg are determined as follows: 


— rt 

* urnttyg) 

V 

■ in' It 

v th 

* IlfV 


By evaluating C lc 


C lc * ’'tv “ C 2 ’'th 


we find that, in general Cj c ^ Cj and therefore, does not lie on the constraint. 

Thus we most vary tg by some value such that C lc goes to Cj . In order to de- 
termine the sensitivity of Cj c to (Ganges in t. we repeat the call to the Conic 
Retired Velocity Routine with tg perturbed Sr a small 6 1 . In a manner similar 
to before a new point (v*^ , v*^ ) is found on Figure 4-3. The factor 


cr = 


C 1 - C »c 

<T *t. - V - c 2 •’“a - V 
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represents the required change in C lc (necessary to extrapolate and v^ to 
die constraint line) divided by the change in experienced by perturbing die 

entry interface time by fit . It is used to extraimlate values for the entry interface 
time, t^ and desired cutoff velocity, v^, which should result in near satis- 
faction of die constraint. 


Wv) ' ^(previous) +<r ^ 

— d(new) ~ — d(previous) + ff fd “ — d (previous) ^ 


The accuracy of this technique depends upon hour much the region over 

which die extrapolation takes place varies from the linearised assumptions. Each 

guidance cycle, however, should bring the extrapolated point closer to the constraint. 


4.8.4 Revising v^ 

Changes in desired terminal velocity, v ^ , affect the predicted terminal 
position (and desired terminal velocity) on die subsequent gtn lance cycle. If this 
effect is ignored, die desired terminal velocity is over corrected and a small 
oscillation in desired terminal velocity is induced. To elii/iinale this over correc- 
tion, a damping factor, p , is introduced for modes 4, 5, 8 as follows: 




* * V„ ) 

—go —go 


where v* 

-go 


is the new velocity-to-be-gained assuming no damping and is given bv 



- ~ — gr&v + —bias 


The damping factor is computed by making use of die approximate partial 
derivative of desired velocity with respect to terminal position 





0 


0 



0 

0 




where t^ 
that 


t -t is the coast time from cutoff to entry interface, 
go 


0 = 


1+ 


-go 


* <*t - * - V 


It can be shown 
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Note that initiall y 0 < p < 1 and p goes to one at the end of the born. 


by 


Velocity-to-be-gained for the neat pass through the routine is now revised 


v _ « ▼* +A»„ 

—go —go —go 


4.9 Throttle Commands 

The throttle setting. . for the Space Shuttle Main Engine and the engine- 
on/ off commands are issued by the last and as yet to be determined block of 
equations. The throttle setting will be assumed to be input each cycle for the ac- 
celeration limited phases or will be computed in die req ui red velocity block of the 
previous cycle for Mode 2 (Reference Trajectory). The array of switches, s «jsme i * 
Sqms i * an d ®nCS i * be utilised to determine which engines will be operating 
i n each phase. 

During the prethrust call, on the last iteration to convergence engine igni- 
tion commands for ^ will be issued. The appropriate time to issue these com- 
mands will be determined by die conditions: s __ = 0, At - 0. 

pre 
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5. DETAILED FLOW DIAGRAMS 

This section contains detailed flow diagrams of die computations and logic 
used for Unified Powered Flight Guidance. The overall flow is illustrated on the 
first diagram (Figure 5- Dwhich shows the nine blocks of computations connected 
by die necessary logic to form a single routine. UPFG. One entry and one exit 
point exist for UPFG with internal branching, as indicated, to perform prethrust 
or active guidance for the various maneuver modes. 

The following diagrams (Figures 5-2 through 5-10) further detail each of die 
uine blocks of in-line computations of UPFG. 
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Figure 5- la. Main Routine 
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ENTER Block 1 



EXIT Block 1 

Figure 5-2. block 1 - Initialization 
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ENTER Block 2 



Figure 5-3a. Block 2 - Update 
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ESTER Block 3 















ENTEH Block 4 



EXIT Block ■» 


Figure 5-5. Block 4 - Integrals of Thrust 
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ENTER Block 5 



Figure 5-6. Block 5 - Turning Rate 


if X - 0 . 
unit ( X > mav cause 
•*o:nputationnl dif- 
ficulties. 


5 - 




ENTER Block 6 


Determine unit thrust Erection b^>ed upon: 
i f * unit(X - iJ/L)>] 



Em* Block 6 


Figure 5-7. Block 6 - Steering Command 
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ENTER Block 7 



EXIT Block 7 


Figure 5-8. Block 7 - Gravity Effects 
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ENTER Blocks 



Figure 5-9a. Block 8 - Velocity- To- Be -Gained 
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EXIT Block 8 

Figure 5 -9c. Block 8 - Velocity-To-Be-Gained (ConU ) 
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Determine engine-off time as function of 
(t + t ) and engine characteristics (TDD) 

s engoff ” 1 



Determine throttle setting for "g - limiting'* (TRDl 

* i 

1 Issue throttle command 



Issue engine-on command, t. 
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EXIT lilui k !) 


Figure 5-10. Block 9 - Throttle Commands 
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6. SUPPLEMENTARY INFORMATION 

Several details of the UPPG program require further study or definition. 

They are listed below, not necessarily in order of importance, 

6. 1 Targeting Assistance 

Since the deorbit maneuver may be a long, low acceleration burn with a 
single OMS engine, traditional targeting techniques based upon an impulsive ma- 
neuver are not adequate. The powered flight guidance routine must assist the 
targeting to determine the Av required for the finite thrust maneuver and compute 
the optimal ignition time. Thus an additional logical path through ITPFG program, 
very similar to the prethrust path, should be included in future revisions. This 
path could be used iteratively by the deorbit targeting program to search for the 
optimal ignition time resulting in minimum Av, 

6.2 Retura-to-Launch-Site Abort 

Considerable effort will be required to define the thrust cutoff condition for 
successful RTLS abort maneuvers. Requirements on dynamic pressure at ex- 
ternal tank separation, orbiter glide back capability, maximum loads, maximum 
vehicle turning rates, fuel depletion requirements, and other factors enter into 
the development of RTLS powered maneuver constraints, 

6. 3 Ascent Initialization 

The problems of efficiently initializing UFFG during an ascent maneuver 
have not been addressed in sufficient detail. It is probably desirable to include 
the boost phase guidance in the overall UPFG scheme. This should simplify 
somewhat the initialization and transition from the relatively simple atmospheric 
phase guidance to the explicit LTG concept used in UPFG, The problems associated 
with changing maneuver objectives during ascent due to engine failures have not 
been explored. The best method of switching from a standard ascent maneuver to 
some abort mode must be determined. The reader should note that the prethrust 
process described in this report is primarily tailored for the on-orbit maneuvers. 

6. 4 External Tank Disposal 

The implications and requirements of the external tank disposal require 
further study. In this document, it is assumed that tank separation will take place 
during a short coast phase after the ascent velocity-to-be-gained has been driven 
to a prespecified value ( < 150 fps). Then the OMS engines complete the ascent 
maneuver. If this technique does not insure a sufficiently accurate external tank 
impact, then further refinement will have to be considered. 
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If the OMS maneuver following tank separation is sufficiently large ( * 300 fps) 
such lhat the length of the OMS maneuver is comparable in length to the phase from 
SRM separation to tank separation, then tne guidance equations must be modified. 
The multi-phase LTG equations, which result in a linear vehicle pitch rate, are 
not optimal for an acceleration profile with a high initial acceleration and a very 
low final acceleration* The final OMS maneuver could probably be treated as a 
constant attitude phase without unduly complicating the equations. 

6. 5 Compensation for Non-Keplerian Gravity Effects 

No attempt has been made to include a technique for compensation of non- 
Keplerian gravity effects in this initial version of UPrG. Further work is re- 
quired to determine whether this compensation should be accomplished during 
premaneuver targeting or during the maneuver. A combination of premaneuver 
compensation with small adjustments during the maneuver will probably produce 
good accuracy with i .limal code* 

6.6 Ascent Maneuver Phase Changes 

The ascent maneuver is divided into several maneuver phases, which may 
include a constant thrust phase, a constant acceleration phase, a final OMS phase, 
and possibly a phase with an assumed SSME failure. Transition from one phase 
to the next may be a function of acceleration, time, or velocity- to-be-gained. 

To solve for several integrals of the thrust acceleration over the total maneuver 
time-to-go, the times at which phase changes occur must be prespecified or cal- 
culated. Since mission perturbations may alter these times, it may be desirable 
to modify them during the maneuver based upon sensed acceleration. 

6.7 Steering Commands 

The development of steering equations to generate commands for the control 
system is incomplete. The steering equations will combine desired vehicle atti- 
tude, desired vehicle rate, and sensed acceleration to produce n steering command. 

6. 8 Throttle Commands 

An algorithm to estimate current acceleration aid calculate an engine throttle 
setting for the constant acceleration ( g -limited) phast of ascent has not been 
developed. The frequency with which these commands must be issued, to pr event 
a sawtooth profile, and the engine response to commands should be evaluated to 
develop this algorithm. 

6. 9 Vehicle Mass Estimate 

The n aintenance of the estimated vehicle mass should be accomplished ex- 
ternal to the UPFG program since this information is also required by the control 
system software. In addition, step changes in the vehicle mass due to events 
such as tank disposal or satellite deployment will obviously have to be handled 
elsewhere. However, for completeness, an equation has been included in this 
document (Block 2) to decrement vehicle mass for simulation purposes. 
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6. 10 Engine Failure 

The failure of an engine during any mission phase could possibly be detected 
through changes in sensed acceleration, however it is assumed that more reliable 
information will be available from the performance monitoring system. Therefore, 
the UPFG equations have been written with the assumption that an external routine 
will notify the UPFG routine of changes in engine status: through the input switches. 

6.11 Tlirottle Lag 

The ascent to reference trajectory guidance mode requires use of the SSME 
throttles. To make the proper adjustments in engine throttle settings, the cur- 
rently commanded throttle settings are used in combination with nominal engine 
performance data. To compensate for engine performance perturbation and throttle 
non parity, it may be desirable to use sensed acceleration information. However, 
this will introduce problems due to both guidance computational delays and eng ; ne 
throttle response. Thus this problem requires further study. 

£yV2 Alternate Input Schemes 

The input list for each maneuver mode, described in Section 3 # has been 
designed with the intent to minimise the size of the flight program. A certain 
amount of ground calculation is required for the ascent guidance modes* however 
on-orbit input is consistant with current CSDL targeting concepts* Future revisions 
could include alternate input schemes for ascent maneuvers which may be easier 
to use but place additional calculi /.on burdens upon the guidance computer. 
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